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O ,■ The purpose of this paper is both to provide mathematical reinforcements to the 



paper [Mecozzi and Bellini : arXiv:1110.1253 [hep-ph]] by taking decoherence into 
consideration and to present some important problems related. 
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_j- ■ We claim that neutrinos have superluminality as a latent possibility. 
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1 Introduction 

k> ■ In September 2011 we encountered a remarkable and unbelievable paper by the OPERA 
collaboration [1] that the speed of neutrino exceeds that of light in vacuum. They 
measured a collective speed of mu-neutrino flying from CERN to Gran Sasso Laboratory ( 
see for example the Fig. 5 in pQ). However, this result conflicts with special relativity in the 
most basic sense. 
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After the paper appeared in the arXiv researchers in the world published or are preparing 
many papers on this topic. Some of them agree with the result, while others don't agree 
according to each author's conviction of the special relativity. See the hep-ph in the arXiv. 

We cannot conclude whether this experiment (detection) is correct or not at the present 
time. Of course, it must be checked by other experiment teams. 

Some researchers protested against the OPERA detection strongly. One of main reasons 
is due to the Kamiokande detection in 1987 [2]. It detected lights and neutrinos coming from 
the Supernova SN1987A at almost the same time. If the speed of neutrino is faster than 
that of light in vacuum, it must have detected neutrinos several years earlier. 

By the way, Mecozzi and Bellini in [3] 1 ! gave a smart interpretation of the result. They 
suggested that the result is due to the superluminal group velocity of neutrinos arising 
from superposition (namely, the neutrino mixing) in Quantum Mechanics [I]. The neutrino 
mixing which is well-known in particle physics is just a quantum mechanical phenomenon. 

However, coherence in Quantum Mechanics is affected by environments and it could 
be destroyed in short time. A long-distance flight from SN1987A might have destroyed 
coherence of neutrinos, and as a result the superluminal group velocity was lost. The paper [3] 
offers an interpretation that the OPERA result does not conflict with Kamiokande detection. 

In this note we provide mathematical reinforcements to the paper [3] in terms of de- 
coherence and would like to offer a "super-smart" interpretation to the OPERA result 

2 Superluminal Group Velocity of Neutrinos 

In this section we review the paper [3] in detail (because it is a bit unclear from the mathe- 
matical point of view). 

First, we prepare some notation for convenience. Since we treat a two level system in 
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the following the target space is C 2 = Vectc(|0i), | ^2)) with bases 





Then Pauli matrices {a x , a y , a z } with the identity I2 



I , <J y 



1 \ / -i \ / 1 \ / 10 

, Oy = , O z = , 1 2 = 

10/ \ i / \ -1 / \01 



act on the space. 

The three generations of leptons are 



e \ / /i \ / r 



^ P / \ Vp 1 \ v T 



However, each neutrino is not single but slightly mixed like 

v' = cos Ufj, — sin z/ r , v' T = sin v^ + cos v T . 

This is called the mixing angle in vacuum, which is small enough. Therefore real genera- 
tions are for example 

/1 



4 / V vL 



Note that the mixing matrix 

(cos sin 
(1) 
— sin cos 

is called the Pontecorvo-Maki-Nakagawa-Sakata matrix ([6], [7], [8]). 

Let us start with a model [9], [TU]. Since we treat two neutrinos {y^ u T } in the paper, the 
Dirac equation for the two neutrinos can be reduced to a Schrodinger form written in terms of 
a two component vector of positive energy probability amplitude in the ultra-relativistic 
limit. 

Then the two neutrino system can be mapped to a two-level quantum system with 
distinct energy eigenvalues along with the assumption of equal fixed momenta [9]. 



We set that rrii (z — 1,2) are the neutrino masses and Ei are the energies given by the 
approximation 



2 S 2 4 

mfc mfc 



Ei = \J(cp) 2 + (mid 2 ) 2 ~ cp-\ — ^ — = cp + 



2p 2pc 

(pi = P2 = p from the assumption) and is the mixing angle. The Hamiltonian becomes 

H = H(p) = (cp + |) 1 2 + e -{sin{2Q)a x - cos(20)a 2 } (2) 

where 

{m\ + ml)c 4 (m 2 , — m\)c 4 e 

£0 = 2p~c ' £ = We ' E ° = C P + T 

Here, note that E — | = E\ and E + | = E 2 . 

Let us rewrite H in (|2]) in a familiar form by making use of Pauli matrices above 

H= . cp + f - § cos(26) fsin(26) 

fsin(26) cp + f + § cos(29) 

For this it is easy to obtain the eigenvalues 

E ± = cp+ 6 j± e - = E ± e - (4) 

and the corresponding eigenvectors 

COS0 

|0_) = | =cos6|0 1 )-sine|0 2 ), 

— sin0 

sinO 

|0+) = | =sine|0 1 ) + cose|0 2 ). 

cos 6 

If we define a unitary matrix 

(cos 9 sin 6 
I = R(Q) 
— sin 9 cos G 

then we can make H diagonal like 

H ---- R(e) [ ) i?(6) f = R(Q) | * ] R(ef (5) 

E + § / \ E 2 



or in a spectral decomposition form 

H=(Eo- I) |0_)(0_| + (Eo + |) |0+)(0+|. (6) 

These forms are used in the next section. 

Next, we label the Hamiltonian by the momentum p. Namely, 

U = H(p)=H(p)®\p)(p\. (7) 

This is in a certain sense the graph of a function. 

If we define eigenvectors of H in ([7]) as simultaneous ones of both flavor and momentum 

H®l\<j)±,p) =E±\<j>±,p), l®p\<p ±1 p) =p\</>±,p), 

then we have 

cos 6 \ 

\<j>-,p) ® \p) = cos9|0i,p) -sin6|02,p), 

-sinG J 

sin0 \ 
\(/>+,p) <8> \p) =sin®\(f) 1 ,p) + cosQ\(f)2,p). (8) 

cos I 
Therefore, the spectral decomposition of H is given by 

H = (Eo - |) \<f>. )P )(<f>.,p\ + (E + |) \<j> + ,p)(<f> + ,p\- (9) 

Next, we consider a time-evolution of the system. 

|-0o> \A) 

I 1 >■ time 

t =0 *=* 

CERN LNGS 

The state \ipt) at time t is given by 

hk)= [dp'e- im ^ h |Vo> (10) 



with the initial state \ipo), and straightforward calculation with (Q gives 

|&> = [ dp'e- UE °/ h {e-^i^p'l^l^p') + e u ^(<P-,p'\M^,p')} (11) 



where E = E (p'), e = e(p'), eo = eo(p') for simplicity. 
From now on we assume some conditions : 

(a) we start with one neutrino flavor, namely, (02 ) p|'0o) — 0, 

(b) the initial amplitude of the neutrino waveform is (0i, p\^o) = (p\ipo)- 
These assumptions seem to be natural. 

Then it is easy to see that (II ip can be rewritten as 

m = f dp'e- itEo ' h {p'\^) x 

[e~ u ^ {sm 2 Q\(j) U p f } + sinG cos Q\<f> 2 ,p')} + e u ^ {cos 2 Q\<f>i,p') - sin0 cos&\4>2,p'}}] 

(12) 

by use of ©. 

Now let us start detection of neutrino : 

(a) first, we perform a flavor measurement (for example, flavor 1), 

(b) immediately after (a), we perform a position measurement. 
When flavor 1 is detected the collapsed state becomes 

\A) — > Wt) = ^jdp'\<t )l ,p')(<f )l ,p'^ t ) = -j=P\ilH) (13) 



where T> is the normalization factor given by 

v = j djf^ t \<j> x ,jf){<j> u p'\ik) = mv\^t) (14) 

«VM) = 1) and 



V = J dj/lfaflifaj/l 

is the projection operator to the flavor 1 state (: V 2 = V, V^ — V). 

It is in general difficult to perform a position measurement immediately after flavor 1 is 
detected, so we average the positions of neutrinos. We believe that this replacement is not 
so bad. 
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The expectation value of position measurement on the collapsed state \ip' t ) is given by 



This is a kind of definition. Note that 



is a position operator, because 



op 



in£-h\x) = ih-^-e- ixp/H = xe~ ixp/h = x(p\x). 
op op 

We must evaluate the expectation value of position (x t ). From f)13p and noting the 
formula 

(<f>i,p'\<f>t : p") =5(p'-p") (15) 

we have and set 

i r , f.^d\ m 



(xt) = pj dp'^tl^p'} [ih—\ (0 ljP U) = p (16) 

for simplicity. 

First, let us calculate (<t>i,p'\i/>t)- From (Tl2|) and the formula (fT5j) it is easy to see 

(<Pup'm = (p>o)F(p') (17) 

where 

F{p') = e - itEa ' h (sin 2 Oe- ite/2h + cos 2 Qe ite/2h ) . (18) 

Then the normalization factor T> in (Tl4|) becomes 

v= [ dp>\(p'\M 2 \F(p')\ 2 - 



Now, we make another assumption. The initial distribution of the neutrino momentum 
Kp'l^o)] 2 is narrow with respect to F(p'), and centered on p' = p (p is fixed). Namely, 

\{p'\iPo)\ 2 *6(p'-p). (19) 

We believe this one natural. Then 

V^\F(p)\ 2 . (20) 
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For later convenience we calculate (I20I) . Here is an elementary formula 
Formula for a, (3 G R 

\ae~ ie + (3e ie \ 2 = (a + (5) 2 - Aaf3 sin 2 9. 



This gives 



\F(p)\ 2 = \sm 2 Qe~ lte/2h + cos 2 Qe ite/2h \ 2 

= (sin 2 6 + cos 2 6) 2 - 4 sin 2 6 cos 2 G sin 2 ( || J 

= 1 — (2 sin 6 cos 6) sin^ 



>2 • 2 ' ^ e 



2h 



= l-sin 2 (2e)sin 2 ^ N ) (21) 

by use of (fTgj) (note that e = e(p)). 

By inserting the equation (|T7|) into A/" in (|T6|) we have and set 

M = Jdp'tyt^p') (ih-^) {^p'^t) 

dp'F{p')(Mp') (ihjj£\ (WoW)} 

dp'\F{p')\ 2 mv') (ih-^j (p'lVo) + J dp'\(p'\M 2 F(p') U-^\ F{p>) 
= Afi+Af 2 

for simplicity. Next, let us calculate A/i and A/2 separately 
From 

A/i = ReA/i + ilmM 



we have 



ImA/i = I(M-M) 



^ J d P '\F(p')\ 2 j(ifo|j/)^Wb) + Wo^W) 

^ J dp'\F(p>)\^\W )\* 
(integration by parts) 
h f dp'\{p'\^)\ 2 ^\F{p')\ 2 

\fdp'5{p'-p)^F{ P r = -- 2 l p \F{p)? (22) 



by use of the assumption in (1191) . Similarly, we have 
RaA/i = ^ (M + A/1) 

= ~ y ^WpOI 2 {<V>ob'> U^\ (p'\iPo) - Wo) (ih£\ (Hp')]- 

The range of integration is narrow enough because of the assumption in (I19p . so we approx- 
imate the integration like 



RaA/i « |F(P)| 2 \fdp'\^ Q \p') (if^j (AM - (p'l^o) (it^ ) (Mp') 
(integration by parts) 

= \F(p)\ 2 J dp' ^o\ P ) U-^j {p'l^o) 



(xo). (23) 

Therefore, we obtain the approximate value 



M = RaA/i H-ilmM « |F(p)| 2 (x ) - vI"I^(p)| 2 - ( 24 ) 

2 ap 



For A/o we have 



A/- 2 = J dp'\(p'\^o)\ 2 F(p') (^ih-^j F(p') 



F(p) [ ih^-\ F(p) = l hF{p)^-F{p) (25) 
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by use of the assumption in ffT9|) . 
Then fl2U) and Q25D give 

= l^)| 2 ^o) + y{Fb)^(p)-^b)F(p)} (26) 

and ( TT6|) and (120]) give (the approximate value) 

(xt) = p = w = (xo) + T ^W ' (27) 

Next, let us calculate the right hand side of fl27|) by use of ( TT8]) : 

F(p) = e" ifi ^ (sin 2 9e- iie / 2ri + cos 2 Qe u ^ 2h ) . 

Noting E = E {p), e = e(p) we have 

|_ F ( p ) = _4^F(p) - i^e- itE °/ h (sin 2 6e- fe / 2S - cos 2 6e^ 2fi ) 
ap n op 2nop x ' 



and 






_ z ^|£ ( sin 2 @e -^/2fi _ cos 2 Qe tie/2H^ ^2 Q( ,Ue/2n + ^2 Qg-ft^R) 



2hdp 

= Awwr-'m^ 9 -'*' 2e+ "» 

where ** is the terms which will be neglected at the final stage. This gives 

ih ( -, s d „. . \ S-En , „, mo t de , 9 ^ 9 ^x 

"2"( F(P) ^ F(P)_C - C J = 'V I(P)I + 2^ (sm 6 ~ C0S 6) 






and (xt) in (1271) is given by 



, . , . 8E t cos(29) de 
10 



Definition The group velocity v g is given by 

/ \ / \ . , \ x t) ~ \ x o) 

Therefore it becomes 



dEp 1 cos(29) de 

Vg dp 2 |F( P )| 2 dp [ ' 

from the result above. 
Remembering 

e (m? + m\)c A (m 2 , — m\)c A 
Eo = CP+ 2' £0 = 2^^' £= 2^ 

from ([2]) simple calculation gives 

dE e <9e e 

<9p 2p y dp p 

and by inserting the above into (129|) we have 

u, = c - — + S— (30) 

9 2p 2p y ' 



cos(26) cos(26) 

|F(p)|2 "" l-sin 2 (2e)sin 2 (|) l J 



where 5 is given by 

from (12TD. 

As a result we obtain 

Theorem (Mecozzi and Bellini) The group velocity v g is given by 



2p + l-sin 2 (2e)sin 2 (||)2p- {6Z) 



Note that the term S was obtained from a quantum effect (the neutrino mixing) and this 
plays a definite role in Superluminal Group Velocity. 
Let us analyze the theorem. For the purpose we set 

_ • *(te 

a = sin — - 

\2h 
11 



to look for some condition satisfying S > 1. Namely, 



s _ cos(26) 

l-asin 2 (26) 
cos(26) > 1 - a sin 2 (26) = 1 - a(l - cos 2 (20)) 

(1 - cos 2 (26))a > 1 - cos(29) (excepting 1 = cos(26)) 

1 1 

a > ; ^ > -• 

1 + cos(26) 2 



Therefore, we obtain 



(|)>i^5>l. (33) 



Here, we assume m 2 3> m\. Then e ~ 6q from ([2]) and S > 1 gives 

v fl wc+(S-l)^>c. 

As a result 



Corollary 1 Under the conditions a > \ and m 2 *^> mi we have the superluminal group 
velocity 

v g > c. (34) 

Note In general, the higher the generation, the heavier corresponding mass. Therefore, the 
assumption m 2 ^> m\ is not unnatural. 

If O = (no neutrino mixing) then S = 1 from (I3"T1) and we have 

eo eo 

under m 2 ^> m\. As a result 

Corollary 2 If = we have the usual velocity 

v g = c (35) 

under a > \ and m 2 >mi. 
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We cannot help admitting a mechanism which accelerates neutrinos arising from the 
neutrino oscillation. How do we interpret the result ? What is the relation to special 
relativity ? In the last part of the paper [3] they write : 

"Of course, this does not mean that the speed of a possible signal transmitted with a neutrino 
wave-packet exceeds the speed of light, it is just a property that comes from the wave-packet 
deformation caused by the interference of the two possible quantum paths that a neutrino 
may follow before reaching the detector" . 

Unfortunately, the author cannot understand what they meant and therefore present the 
following 

Problem Give a mathematical expression to their claim. 

3 Decoherence of the Neutrino Oscillation 

In this section, for the model in the previous section we take decoherence into consideration 
in order to make it more realistic. Then we can build a bridge between Corollary 1 and 
Corollary 2. For a general introduction to this topic see for example [TT] . 
First of all we present the following 

Problem Is there no problem to apply theory of decoherence to neutrinos in a long-distance 
flight ? 

Although this problem is very subtle, let us proceed to the discussion of decoherence. 

Since the two neutrino system can be mapped to the two level system we prepare some 
notation from Quantum Optics. For 

1 \ ( 



it is easy to see 



0/ Z V 1 



1 \ / 



0/ V01 
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Let us remember 



Ei \ T 

H R(G) ( R{ef 

E 2 ) 

E x cos 2 6 + E 2 sin 2 6 (E 2 - E 1 ) sin Q cos 6 
(E 2 -E l ) sin6 cos6 £ x sin 2 6 + £ 2 cos 2 Q 



(36) 



from (E]) and set 



Ex 
flb = | | (37) 

E 2 



Note that H and f/o are symmetric matrices (H = H T , H = Hq). 

To treat decoherence in a correct manner we must change models based on from a pure 
state to a density matrix. The general definition of density matrix p is given by both p* = p 
and trp = 1, so we can write p = p(t) as 

a b 
p=\ ) (a = a, d = d, a + d = 1). (3£ 

b d 
The general form of master equation is well-known to be 

— p = —i[H, p] + Dp (■<= /I = 1 for simplicity) (39) 

where 

( 1 \ \ ( 1 1 

Dp = n I ^-P^+ - 2 ff + a -P - 2^°"+°'- J + v [ a +P a - ~ 2 a - a +P ~ 2p a - a + 



I- 



and /i > v > q We must solve the equation. 
If we write H in ([36]) as 



h k 
H= | (h, k, I e R) (40) 

fc / 



2 we don't know how to determine the precise values of [i and v in our system 
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for simplicity, then the master equation above can be rewritten as 





fa ) 


( 


d 


b 




dt 


b 






K d ) 


\ 



—ji ik 

ik i(l-h)-^ 

—ik 



-ik 



—ik 

-i(l - h) - 

ik 



2 



„ \ 




fa\ 


—ik 




b 


ik 




b 


-" ) 




IV 



(41) 



We leave the derivation to readers. 
Note and set 

—[/, ik 

ik i{l-h)-^f 



-ik 




v 
-ik 



-ik 

fj, —ik ik 

ik -ik 

ik i(l — h) —ik 

—ik —i(l — h) ik 

-ik ik 



-i{l -h)-^- ik 



-v 



\ 



The general solution of (14T1) is given by 

' a(t) ^ 
b(t) 
b(t) 

\ d(t) J 



( 



+ 



-/'■ 








V 





fl+U 

2 














fl+U 

2 





u 








— 1 



\ 



V 



H + D. 



JIH+D) 



' tt(0) ^ 

6(0) 

5(0) 

V d(o) y 



(42) 



However, it is not easy to calculate the term e*^ + > exactly, so we use a simple approxima- 
tion 

e t(H+D) _ e t(D+H) _ e tD e tH_ 

In general, we must use the Zassenhaus formula, see for example [12], [15] . Therefore, we 
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treat the approximate solution 



e tD e tH 



/ a(t) ^ 
b(t) 

b(t) 

\ d(t) J 



First, we calculate e tD . For the purpose we set 



1 o(0) ^ 
6(0) 

5(0) 

v d(0) / 



(43) 



K 



—fj, v 
/I — v 

and calculate e tK . The eigenvalues of K are {0, — (/i + z/)} and corresponding eigenvectors 
not normalized) are 

'A ( ' 

, -0 + i/) < — ► I 



If we define the matrix 



O 



then it is easy to see 



v 1 

H —1 

K = 



o- 1 



(JL + V 



1 1 



/i — ^ 



o- 1 



■{» + ") 



and 



J-ftT 



O 



o- 1 



v + fie-*^-^ v - i/e - *^") 



Therefore, we have 



„t£> 








(j,— fie 



-t(n+i>) 



\ A*+" 



-t(/i+;/) 





e 2 






f x + u \ [1- /ie-'^ +ly ) /i + ue~ 



-Uli+v) 





e 2 








M+lye -t(M+^) 



/i + ^ 



z/ v » 





. /i /i , 



if t is large enough (t ~^> 1/V). 



(44) 
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Next, we calculate e tH . Since we need some properties of tensor product in the following 
see for example [13]. We can write the equation as 



H = -i (H <g> 1 2 - 1 2 ® H) (^H = H 7 



In fact, 
H = 



h k \ / 1 

® 

fc i / \ 1 

r / 



-« < 



It is well-known that 



h 





k 





/ 





h 





fc 




k 





I 










k 





'J 


\ 





. 


1 







1 





1 


ft 


A- 








A- 


/ 














h 


k 








k 


I 



® 



h k 
k I 



\ 


/ 




> = -i 


J 


\ 



A- 




-A; A; 

■(l-h) fc 

l-h -k 

k 0-fc 



,tH __ e -it(if(g)l2-l2®H) __ e -itH®l 2 & itl2®H _ ( e -itH gj Wj ^ gitff\ 



e -ftH g, gttff 



Since 



H = R(Q) | * | R(ef 

£ 2 



we have 



e~ itH = R(Q) 



-itE-i 



e 



-it£ 2 



«(e) 5 



and 



-itEi 



AH 



R(Q) 







e~ itE2 

(R(e) <g> i?(e)) 



i2(0) 



T 



; iiJ5l g 



i?(0) I I fl(6) r 

e if£2 



it£i 



e" iijB2 



/ 



{R(G)®R(Q)) 



e itE 1 Q 

e iijB2 
\ 



(R(Q) ® i2(0)) J 



e it(£ 2 -Bi) 



-it(E 2 -E 1 ) 



(R(B) ® R{Q)f . 



1 
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Since 



R(G) 



cos B sin 
— sin B cos B 



we have 



/ 



R(G) <g> R(G) 



cos B sin B cos B sin B 
cos 2 B - sin 2 B 

- sin 2 B cos 2 B 



sin 2 B 
cos B sin B 
cos B sin 6 



cos 2 B 

— cos B sin B 

— cos B sin B 
sin 2 B — cos B sin B — cos B sin B cos 2 B 

and, by setting J = e lt ( E2 ~ El ^ = e tte for simplicity, 



\ 



where 



and 



JH 



R(Q)®R(Q) 



J 



J 



\ 



(R(G)®R(G)f 



1 



( \ 

«11 «12 «13 «14 



V 



«41 «42 «43 «44 
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a u = cos 4 6 + (J + J -1 ) cos 2 9 sin 2 9 + sin 4 9, 

ot\2 = — cos 3 B sinB + Jcos 3 BsinB — J -1 cosBsin 3 B + cosBsin 3 B, 

«i3 = — cos 3 B sin B — J cos B sin 3 B + J -1 cos 3 B sin B + cos B sin 3 B, 

«14 = cos 2 B sin 2 B - ( J + J~ x ) cos 2 B sin 2 B + cos 2 Bsin 2 B 

a 4 i = cos 2 B sin 2 B - ( J + J" 1 ) cos 2 B sin 2 B + cos 2 Bsin 2 B, 

a 42 = — cos 6 sin 3 6 — J cos 3 B sin B + J~ l cos B sin 3 B + cos 3 B sin B, 

a 43 = — cos B sin 3 B + J cos B sin 3 B — J -1 cos 3 B sin B + cos 3 B sin B, 

a 44 = sin 4 B + (J + J" 1 ) cos 2 Bsin 2 B + cos 4 6. 
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Note that *'s in the matrix are elements not used in the latter. We leave this derivation to 

readers. 

Here, we list very important relations among {a} 



an + a 4i = 1, ai2 + a 42 = 0, a i3 + a 43 = 0, a u + a 44 = 1. 



Therefore, from (T43]) and ()45p we obtain 

/ . n n .. 



1 a(t) ^ 
b(t) 

b(t) 
d(t) f 



fi + v 



fJ> + v 



v v 





y fi ft J 

v v 





. fi /! , 



an ai2 ai3 au 

* * * 

* * * 



* 



* 



a 4 i a 4 2 a 4 3 a 44 
\ 



( o(0) ^ 
6(0) 

5(0) 



o(0) 
K0) 
6(0) 

V d(o) y 



for t > 1/za 



For H in (J3ZJ) 

' Ei 

E 2 
we can perform the same process much easily. The master equation is 



H f) 



( \ ( 

a 



d_ 
dl 



-H 

i(l-h)-^ 


V t* o 

and the (exact) solution is given by 



\ d / 



-0 v 



-i(l -h)-^ 

-v 



\ 


fa\ 




b 




6 


J 


w 



(45) 



(46) 



1 a(t) ^ / 



6(t) 

6(t) 

v d(t) y 



fl+U 






fi—fie 
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When i> \jv we obtain 
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As a result we have 

Theorem Two systems (master equations) whose Hamiltonians are H and H have the 
same asymptotic behavior ( l4"6j) and ( 14"T|) under our approximation. 

This theorem implies 

Corollary The mixing angle G will become if t is large enough. 

From both this corollary and corollary 2 in the preceding section we can conclude that 
the speed of neutrinos (after a long-distance flight) is just that of light in vacuum. 

By the way, our calculation in this section is based on a simple approximation. This is a 
bit poor, so we present the following 

Problem Give the explicit (full) calculation. 

As for interesting topics of decoherence (which is essential in Quantum Physics) arising 
from Quantum Optics or Quantum Computation see our papers [H], [15] and [16], [T7| . 

At the end of this section, one comment is in order. It seems to the author that Neutrino 
Physics gets along with Quantum Optics or Quantum Computation, see for example [T8] . 
In order to make some (deep) relations clear further studies will be required. 

4 Concluding Remarks 

In this paper we re-examined the paper [3] by Mecozzi and Bellini in detail and tried to give 
mathematical reinforcements to it by taking decoherence into consideration. Our conclusion 
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IS 

Neutrinos are latently superluminal. 

We would like to present the following 

Problem Re-check our result from a different point of view. 

Let us write once more that our argument is based on group velocity. Therefore, it is 
a bit unsatisfactory. 

Whether the OPERA experiment is correct or not is not concluded at the present time 
and it must be checked by other experiment teams. However, such a check will take time. 
Therefore, it is very important for us to state 

Problem Make some (inside) questions clear from a theoretical point of view. 

Regarding papers related to this topic see for example [19] . 

The work is a great challenge to not only (young) Physicists but also (young) Mathe- 
maticians, so we conclude the paper by citing famous sentences by late Steve Jobs u 

Stay hungry, Stay foolish 

(see the Concluding Remarks in [S]). 

Acknowledgments We wish to thank Shin'ichi Nojiri and Ryu Sasaki for useful suggestions 
and comments. 
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